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We present a derivation of Boltzmann principle Sb — fcalnW based on classical mechanical 
models of thermodynamics. The argument is based on the heat theorem and can be traced back to 
the second half of the nineteenth century with the works of Helmholtz and Boltzmann. Despite its 
simplicity, this argument has remained almost unknown. We present it in a modern, self-contained 
and accessible form. The approach constitutes an important link between classical mechanics and 
statistical mechanics. 
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I. INTRODUCTION 

One of the most intriguing equations of modern physics 
is Boltzmann's celebrated principle 



S B = fcslnW, 



(1) 



where ks is Boltzmann's constant. Despite its unques- 
tionable success in providing a means to compute the 
thermodynamic entropy of isolated systems based on 
counting the number of available microscopic states W, 
its theoretical justification remains obscure and vague 
in most statistical mechanics textbooks. In this respect 
Khinchin commented: [l[ "All existing attempts to give 
a general proof of this postulate must be considered as 
an aggregate of logical and mathematical errors superim- 
posed on a general confusion in the definition of the basic 
quantities." This lack of a crystal clear proof of Boltz- 
mann's principle puts physics students, teachers, and in 
indeed all physicists, in the uncomfortable position of 
being forced to accept the formula as a postulate that is 
necessary to link thermodynamic entropy to microscopic 
dynamics. 

Recent studies in the field of history and foundations 
of statistical mechanics Q have drawn the attention to 
the fact that a similar formula, 



S = fc B ln$, 



(2) 



emerges naturally from classical mechanics if (a) the er- 
godic hypothesis is made, (b) the properties that entropy 
should satisfy are appropriately set and the basic quan- 
tities are consistently defined. The quantity <f> is the vol- 
ume in phase space enclosed by a hyper-surface of con- 
stant energy E. 

Equation (J2J is valid for both small and large systems 
and coincides with the Boltzmann formula for large sys- 
tems. Hence the derivation of Eq. © provides the miss- 
ing proof of Eq. |T]). The basic argument underlying the 
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derivation of Eq. @ can be traced back to as early as 
the second half of the nineteenth century in the works of 
Helmholtz and Boltzmann. 0, Q 

The purpose of this article is to provide a widely ac- 
cessible modern presentation of the original argument of 
Helmholtz and Boltzmann 0, 0] and of its recent devel- 
opments that could be used in the classroom. Wc 
derive Boltzmann's principle from classical mechanics 
with one simple guiding principle, viz., the heat theorem 
(see Statement 1, below), and one central assumption, 
namely, the ergodic hypothesis. 

In Sec. [H] wc briefly review basics of thermodynam- 
ics. We give concise formulations of the first and second 
laws of thermodynamics and introduce the heat theorem. 
We then construct a one-dimensional mechanical model 
of thermodynamics in Sec. IIHI according to the work 
of Helmholtz. 0] The concepts of ergodicity and micro- 
canonical probability distribution emerging naturally in 
this model are introduced in Sec. IIVI We then generalize 
the one-dimensional model to more realistic Hamiltonian 
systems of iV-particles in three-dimensions in Sec. [V] At 
this stage the ergodic hypothesis is made and Eq. @ is 
derived. In Sec. IVI1 we point out that the mechanical 
entropy of Eq. @ does not change during quasi-static 
processes in isolated systems, in agreement with the sec- 
ond law of thermodynamics. Non quasi-static processes 
that can lead to an increase of entropy have been treated 
elsewhere. 0, 0] In Sec. IVIII the Boltzmann principle is 
derived. A summary and some remarks concerning the 
validity of the ergodic hypothesis arc given in Sec. IVIIII 

In the text we present the line of reasoning and the 
main results, while proofs and problems are provided in 
the appendices. 



II. CLAUSIUS ENTROPY 

The purpose of this paper is to construct a classical 
mechanical analog of thermodynamic entropy. To this 
end it is necessary to give a clear account of the definition 
of entropy in thermodynamics. We now review the first 
and second laws of thermodynamics in the formulation 
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given by Clausius (see Ref. Q). The latter gives the 
definition of thermodynamic entropy. 



First Law of Thermodynamics 

In its differential form the first law of thermodynamics 
reads: Si 



dE = 5Q + SW, 



(3) 



where dE is the change in internal energy, SQ is the heat 
added to the system and SW is the work done on the 
system during an infinitesimal transformation. The first 
law is the energy conservation law applied to a system in 
which there is an exchange of energy by both work and 
heat. 

Of crucial importance for the understanding of the first 
law is that SQ and SW are inexact differentials, whereas 
dE is exact. The internal energy E is a state variable, 
namely a quantity that characterizes the thermodynamic 
equilibrium state of the system. On the other hand, W 
and Q are quantities that characterize thermodynamic 
energy transfers only and are not properties of the state 
of the system. @ 

Exactness of differentials is best understood in terms 
of their integral along a path in the system's state space. 
The differential is exact if and only if the integral de- 
pends only on the end points of the path. Accordingly, 
the integral of an inexact differential along two different 
paths with same end points may take on different values. 
The interested readers find a summary of the formal def- 
inition of differential forms and their major properties in 
Appendix fA] 



Second Law of Thermodynamics 

The second law of thermodynamics can be conve- 
niently summarized as three statements. 

Statement 1. The differential SQ/T is exact. 

Statement 1, is one of the most important statements 
of thermodynamics: Although SQ is not an exact differ- 
ential, when it is divided by the absolute temperature T, 
an exact differential is obtained. 

This is equivalent to stating that that there exist a 
state function S. such that 



6 -®=dS. 
T 



(4) 



The function S is called the thermodynamic entropy of 
the system. 

Statement 1 can be expressed in an equivalent way also 
in integral form, by stating that the integral of SQ/T 
along a path connecting a state A to a state B in the 
state variables' space, does not depend on the path but 
only on its endpoints A and -B.[ll[ This in turn says that 



there exists a state function S (i.e. the thermodynamic 
entropy), such that 



SQ 
T 



= S(B)~S(A). 



(5) 



From Eq. (|3|) it is SQ = dE — SW. In general, the work 
is performed by changing a certain number of external 
parameters Ai, e.g. volume, magnetic field, electric field. 
Then the work SW is given by — J2i Fid-K, where Ft de- 
note the corresponding conjugate forces, i.e., pressure, 
magnetization, electric polarization, respectively. There- 
fore it is: 



SQ = dE + F t dXi 



(6) 



Without loss of generality, in the following we will restrict 
ourselves to the case of only one external parameter V 
with conjugate force P:fl2| 



8Q = dE + PdV 
In this case Statement 1 can be re-expressed as: 



(7) 



(dE + PdV) jT = exact differential = dS, (8) 

This is referred to in the literature as the heat theorem. 
[Io| The heat theorem can be re-expressed in equivalent 
terms as: there exists a function S(E,V) such that: (l3j 



(9) 



dS_ _ 1 

dE ~ T ! dV ~ T' 

It is worth emphasizing that any inexact differential, 
like for example SQ = dE+PdV, does not enjoy the same 
property: it is impossible, in general, to find a function 
of state Q(E, V) such that dQ/dE = 1 and dQ/dV = P. 

The following two statements, regarding the function 
S complete Clausius's form of the second law: 

Statement 2. For a quasi-static process occurring in 
a thermally isolated system, the entropy change between 
two equilibrium states is zero, 



AS = 0. 



(10) 



Statement 3. For a non quasi-static process occurring 
in a thermally isolated system, the entropy change be- 
tween two equilibrium states is nonnegative, 



AS>0. 



(11) 



A crucial point that must not be overlooked is that 
Statements 2 and 3 pertain to thermally isolated sys- 
tems. This means that the system is not in contact with 
a thermal bath, by means of which one could in princi- 
ple control either its temperature or its energy. Thus the 
processes mentioned therein are processes in which only 
the external parameter V is varied in a controlled way 
and there is no control over the variable E. In State- 
ment 2 the change of the parameter V is so slow that at 
any instant of time the system is almost at equilibrium 
(quasi-static process). In Statement 3, this requirement 
is relaxed. 
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III. ONE-DIMENSIONAL MECHANICAL 
MODELS OF THERMODYNAMICS 

In this section we construct a one-dimensional classical 
mechanical analogue of Clausius thermodynamic entropy. 
This construction dates back to Hclmholtz 0, H, Q3 and 
is based on the heat theorem (J8]). 

Consider a point particle of mass m and coordinate x 
moving in a [/-shaped potential <p(x), as illustrated in 

Fig. m 

In order to allow for the possibility of doing work on the 
particle by means of an external intervention, we assume 
the potential <p to depend on some externally controllable 
parameter V: (p = ip(x;V). As an example one could 
think of a pendulum whose length can be changed at will 
by an experimenter while the pendulum oscillates. In 
this case V would denote the length of the pendulum. 
The Hamiltonian of the system is: 



H(x,p; V) = K(p) + cp(x;V), 



(12) 



where K(j>) = p 2 /2m is the kinetic energy and p is the 
momentum. 

Now that our mechanical system is defined, we have 
to specify its "internal energy" E, "temperature" T, and 
the "force" P conjugate to the external parameter V. 

For the internal energy we simply take the energy E 
given by the Hamiltonian. For a fixed V the particle's 
energy E is a constant of motion. For the sake of sim- 
plicity we chose the gauge of the potential in such a way 
that the minimum of the potential is regardless of the 
value of V. 

Once V and E are specified, the orbit in phase space of 
the particle is fully determined: we say that E and V are 
the system's state variables. 14| The particle moves back 



and forth between the two turning points x±(E, V), and 
draws closed orbits in phase space with a certain period 
t(E,V). See Fig. [TJ 

Now that the state variables are fixed we have to define 
the corresponding temperature and conjugate force. In 
agreement with the common understanding of tempera- 
ture as a measure of the speed of the particles, we take 
the temperature to be proportional to the kinetic energy 
averaged over one period: 



T(E,V) 



k B r(E,V) Jo 



■(E,V) 



dtK(p(t;E,V)) (13) 



Thanks to the the factor 1/fcs, T has the correct dimen- 
sions of a temperature. For the conjugate force we take 
the time average of — ^ (HI, i.e., 



P{E,V) :=- 



1 



t(E, V) J 



T{E ' V) dt dip(x(t;E,V);V) 



dV 



(14) 

In Eqs. qi3|14|> . x(t;E, V) and p(t; E, V) are the so- 
lution of Hamilton's equations of motion with a fixed 
V, and an arbitrary initial condition x$,po such that 
H(x ,po;V) =E. 
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FIG. 1: Point particle in the [/-shaped potential (p(x; V) = 
mV 2 x 2 /2. Top left panel: shape of the potential for a certain 

V — Vi . Top right panel: shape of the potential for a certain 

V = Vi. Bottom left panel, phase space orbit corresponding 
to the potential tp(x;Vi) at energy E\. Bottom right panel, 
phase space orbit corresponding to the potential ip{x\Vi) at 
energy E2- The two quantities E, V , uniquely determine one 
"state", i.e., one closed orbit in the phase space. 



Having identified the mechanical analogues of internal 
energy, external parameter, temperature and conjugate 
force with the quantities E, V, T, P, respectively, we can 
now ask whether a mechanical analogue of entropy, S, 
exists. In order to answer this question we must ask, 
in agreement with Statement 1 as expressed in Eq. ([5]), 
whether there exists a function S(E, V) such that: 



dE y ' ' T(E,Vy dV 



^ {E ' V) = T(E^V)- (15) 



The answer to this question is given by the following 
theorem. 

Theorem 1 (Hclmholtz). A function S(E,V) satisfying 
Eq. M5)) exists and is given by 

rx+{E,v) 

S(E,V) = fcslog2 / dxy/2m{E - ip{x;V)) 

Jx-(E,V) 

(16) 

The proof of the theorem is reported in Appendix [Bl 
see also Ref. ( pp 45-46). 

The quantity \/2m{E — ip(x; V)) represents the mo- 
mentum of the particle when it is located at x. Solving 
the equation E = p 2 /2m + ip(x; V), and taking the posi- 
tive root, we obtain: 



p(x] E, V) = y/2m{E - <p(x\V)) (17) 
Thus the entropy can be rewritten compactly as: 

S(E,V) = k B \og<j)pdx (18) 

The integral § pdx is called the reduced action. [l5| It rep- 
resents the area $ enclosed by the orbit of energy E, and 
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parameter V in phase space: 

S(E,V) = k B \og<S>{E,V) 

where 



dpdx. 



H(x,p;V)<E 



(19) 



(20) 



Theorem 1 says that there exists a mechanical counter- 
part of entropy, and that this is the logarithm of the 
phase space volume enclosed by the curve of constant 
energy H(x,p; V) = E. 

The fact that there exists a function S satisfying Eq. 
(fT5| is a highly iron trivial result which tells us that 
we have a consistent one-dimensional mechanical model 
of thermodynamics. Although this model nicely sug- 
gests the deep connection between classical mechanics 
and thermodynamic entropy, it is definitely too stylized 
to model a real thermodynamic system composed of as 
many as 10 23 particles. It is necessary to generalize 
Helmholtz Theorem to multidimensional systems. 



IV. ERGODICITY AND MICROCANONICAL 
ENSEMBLE 

The main ingredients needed for the generalization of 
our model to more degrees of freedom are ergodicity and 
the microcanonical ensemble. 

We use the one-dimensional example of the previous 
section to introduce these important concepts. Imagine 
we want to calculate the time average of a phase function 
f(x,p) over the orbit specified by E and V: 



(f) t :=- [ T dtf(x(t),p(t)). 



(21) 



For simplicity of notation we have dropped the explicit 
dependence of t, x(t), p(t) and (f)t on E, V. 
Since p = mv = mdx/dt the differential dt is: 



dt = 



dx 
p(x) 



(22) 



With this we obtain: 



, , 2m f x + dx , , NX 

where the factor 2 stems from the fact that the particle 
goes from x_ to x + , in a half period, i.e., r/2. Now 
consider the following integral 



(p 2 /2m + <p(x; V) - E) 



(24) 



where 8 denotes Dirac's delta function. Using the formula 
${f(p)) = Y,i S (P - Pi)/\f'(j>i)\, where the p^s are the 
zeroes of f(p), and / dp8(p — Pi)=l, we get: 

/ dp5 (p 2 /2m + <p(x; V) — E) = 2m/p(x) (25) 



Then Eq. J23J becomes: 
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(f) t = -JdxJdp6 (p 2 /2m + <p(x; V) - E) f(p, x) 

(26) 

where the integration extremes x± need not be specified, 
being implied by the Dirac 6. The period r is given by: 

r dt - 2 r + — 

JO ~ Jx- P(x) 



dx j dpS (p 2 /2m + ip(x; V) - E) (27) 
Hence we arrive at: 

(f)t = [dx [ dpp(x,p;E,V)f(p,x) (28) 



where we have introduced the phase space probability 
density function 



p{x,p;E,V) 



1 



-(E,V) 



5{p 2 /2m + Lp{x;V)-E) (29) 



From Eq. (|29|) it is clear that r is the normalization. 
This p(x,p; E, V) is called the microcanonical distribu- 
tion. Eq. (|2"8")l says that the time average of a phase 
space quantity f(x,p) over one period, is equal to its mi- 
crocanonical average. This property is called ergodicity. 
All one-dimensional systems with a [/-shaped potential 
are ergodic. 



V. MULTI PARTICLE MECHANICAL MODELS 
OF THERMODYNAMICS 

We now can extend the previous treatment to systems 
of A-particles in three dimensions. The Hamiltonian for 
an interacting system of iV-particles of mass m is 

H N (<i,p;V)=K N (p) + <p N (q;V), (30) 

where Kn{p) = Xa=iP? ^ s the kinetic energy, tp^ 
is the potential energy, and the coordinates q = {qi}^ 1 
and their conjugate canonical momenta p = {pi}f=i are 
SA^-dimcnsional vectors. 

In analogy with one-dimensional systems with a U- 
shaped potential, wc define the microcanonical probabil- 
ity distribution as 



p N (q,p;E,V) 



1 



6(E-H N (<i,p;V)), (31) 



tt N (E,V) 
where Qn(E,V) is the normalization: 

Q N (E, V) = J ■ ■ ■ J dqdpS (E - H N (q, p; V)) . (32) 

Continuing the analogy with one-dimensional systems, 
we make the following crucial assumption: 
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Ergodic Hypothesis. For given E and V, the time av- 
erage (f) t of any function /(q, p) is uniquely determined 
and is equal to its microcanonical average (f)^, i.e.: 

(/)* = J - J dqdp p N (q,p;E,V)f(q,p) = (f)^ 

(33) 

In analogy with Eqs. (|13ll4[) we define the temperature 



as: 



and the conjugate force as 



^ ■■-{¥), 



(34) 



(35) 



We ask, in agreement with Statement 1 as expressed in 
Eq. ([9]), whether there exists a function Sn{E,V) such 
that 



dS 



N 



dE 



(E,V) 



1 



T N (E,VY dV 



dSN (E,V)= PN{E > V) 



T N (E,V) 

(36) 

The answer is given by the following theorem: 

Theorem 2 (Helmholtz, Generalized). A function 
Sn(E,V) satisfying Eq. 136\) exists and is given by: 



S N (E,V)=k B log$ N (E,V) 



whe 



*n(E,V) 



dqdp 



(37) 



(38) 



H N (q,p)<E 



The proof, which is based on the equipartition theo- 
rem, is given in Appendix [C] See also Ref . . 

We draw the attention to the fact that, unlike temper- 
ature and conjugate force, Sn is not in the form of the 
time average of some phase function /(q, p) 

Theorem 2 says that ergodic systems constitute ideal 
mechanical models of thermodynamics. One can de- 
fine their state variables by E and V as in thermo- 
dynamics. Moreover, one can define their temperature 
and conjugate force straightforwardly as functions of 
the state variables. Surprisingly, the heat differential 
(dE + PnoIV)/Tn is exact, allowing for a consistent and 
logical definition of entropy Sn- 



VI. ADIABATIC INVARIANCE 

According to Theorem 2, Sn complies with the first 
law of thermodynamics and statement 1 of the second law 
of thermodynamics: Is the construction consistent with 
statements 2 and 3 of the second law of thermodynamics 
as well? 

Let us focus on statement 2. In order for Sn to be 
consistent with this statement it is necessary that if the 



parameter V is varied very slowly in time (much slower 
than any time scale of the system dynamics) from a cer- 
tain Vb = V(to) to a certain Vf = V(tf), the correspond- 
ing change of the entropy Sn is null. Note that, by al- 
lowing for a time dependence of V, the system's Hamil- 
tonian now becomes time dependent, and energy is not 
conserved. Let the system be at t = to, in qo,Po- Under 
the time dependent Hamiltonian 



H N (q,p;V(t)) = K N (p) + <p N (<KV{t)), 



(39) 



it evolves to a new phase space point 
Q/(QO)Po)i P/(qo, Po), where we made explicit the 
dependence on the initial condition of the evolved 
phase space point. Then the energy at time tf is 
E±= H N (qj(qo,Po),Pf(qo,Po);V(tf)). It is known 
)l6l | that, for ergodic systems, the energy reached at the 
end of a very slow protocol depends only on the initial 
energy E = HN{qo,Po]V(t )), and is determined by 
solving the following equation for Ef. 



<f> N {E ,V ) = ®N(E f ,V f ) 



(40) 



That is, the quantity $>n(E, V) does not change in the 
course of time, when V is varied infinitely slowly. This 
property is called in classical mechanics adiabatic in- 
variance. Since it is Sn(E,V) = k B log &n(E, V), and 
<S>n(E, V) is an adiabatic invariant, it is evident that Sn 
is an adiabatic invariant too. Namely, it does not change 
if V is changed very slowly in time. Thus Sn complies 
with Statement 2. For completeness in Appendix |Pl wc 
provide a proof of adiabatic invariance of &n- See also 
Ref. (3 (pp. 27-30). 

It is also possible to prove that, in an averaged sense, 
Sn complies with Statement 3, as well.JH, @ In this case 
one has to consider the average change of entropy, be- 
cause, for fast transformations, the final energy is not 
uniquely determined by the initial energy, and depend- 
ing on the initial conditions, one ends up with different 
final energies, i.e., different final entropies. 



VII. BOLTZMANN PRINCIPLE 

For a system composed of a very large number N of 
particles which interact via short range forces, the phase 
space volume <&n(E) approaches an exponential behavior 
<f> N (E) oc e E . Since = d$> N /dE (sec Eq. (|32])L it is 
$at oc n N (see Ref. [13, p. 148). 

The quantity O/v, defined in Eq. (|32[) . represents the 
measure of the shell of constant energy Hn(h,p;V) = 
E. As such it is proportional the number W of micro 
states compatible with the given energy E. (According to 
scmiclassical theory each micro state occupies a volume 
h 3N of phase space, where h is Planck's constant. [HI By 
introducing an arbitrary energy scale AE, the number 
W is given by W = n N AE/h 3N ). Thus, for very large 
N, 



$at oc ft N oc W, N > 1 



(41) 
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By taking the logarithm, we have 

S N ~ fc B lnW = S B , AT>1 



(42) 



except for an irrelevant constant. Eq. (|4"2")) says that 
for large ergodic systems composed of particles interact- 
ing via short range forces, the differential of Boltzmann 
entropy is equal to the differential SQ/T. Hence it can 
be identified with Clausius entropy. This is a proof of 
Boltzmann principle. 



VIII. CONCLUSIONS 

Given an ergodic system, it is possible to specify its 
thermodynamic state by means of the total energy E 
and the external parameter V. Given the state E, V, we 
can unambiguously define the quantities Tn(E, V) and 
Pn(E,V). Once these are identified with the system 
temperature and conjugate force, one can ask whether, 
as prescribed by the heat theorem, the combination: 

dE + P N dV 



T N 

is an exact differential. Surprisingly the answer is posi- 
tive, meaning that there exist a function Sn (E, V) which 
can be identified with the thermodynamic entropy of the 
system. The generalized Helmholtz theorem says that 
this is given by the logarithm of the volume (E, V) 
of phase space enclosed by the hyper-surface of energy 
H (q, p; V) = E. For macroscopic systems this entropy 
coincides with Boltzmann entropy, thus revealing the ra- 
tionale of Boltzmann principle. 

The entropy in Eq. (|37|) is sometimes referred to in 
the literature as Hertz entropy. [2(| Hertz [2l|, [22| de- 
rived it from the requirement of adiabatic invariance (see 
also Refs. 0, [23|, |24|), whereas we have derived it here 
from the heat theorem. Its fundamental character is also 
recognized in Ref. [25| where its property of being a 
canonical invariant is emphasized and in Ref. [2r| which 
highlights its compliance with the equipartition theorem, 
and the fact that it is a positive and increasing function 
of the energy [H HI- The entropy in Eq. J37]) also ap- 
pears in Gibbs seminal book.[29| However its connection 
with the heat theorem has not previously recognized. 

The most crucial point of the derivation of Boltzmann 
principle is the introduction of the ergodic hypothesis. 
Although this hypothesis is generally believed to be true 
for real macroscopic systems, its mathematical proof is 
a formidable challenge which has been achieved only in 
few special cases. [3(J A proof that a gas of elastically col- 
liding hard spheres is ergodic was announced in 1963 by 
Sinai. [3l| However the full proof was not published and 
the problem is still open (see Ref. [HJ for a more detailed 
discussion). Nonetheless ergodicity of hard spheres sys- 
tems seems plausible as indicated also by recent numeri- 
cal simulations (see Sec. IV of Ref. [H|). The hypothesis 
cannot be true in the case of real crystals. Here the nuclei 



remain close to their lattice site, preventing them from 
sampling the whole energy hypersurface homogeneously 
(this being a necessary condition for ergodicity). 

In regard to these difficulties, it is worth pointing out 
that the present derivation of Boltzmann principle does 
not make use of the fact that the average of any arbitrary 
phase function be equal to its micro canonical average, 
as required by the ergodic hypothesis. It only requires 
that the time average of K and —dcp/dV, be equal to 
their microcanonical averages (see the proof of Theorem 
2 in appendix [C| . Thus the ergodic hypothesis can be 
greatly relaxed by requiring the much less stringent con- 
dition that Tjv and Pjy, Eqs. (|34I35|) . can be calculated 
as microcanonical averages. [Hf In this case, the Clausius 
entropy can still be calculated via the Hertz formula (|37[) . 



Acknowledgements 

We wish thank the Texas Section of the American 
Physical Society for the "Robert S. Hyer Recognition for 
Exceptional Research" presented at its Fall 2008 meeting, 
and Cosimo Gorini for reading the manuscript. We also 
would like to thank Prof. Randall B. Shirts for providing 
the translation of Ref. [25| . Valuable remarks from the 
anonymous referees are gratefully acknowledged. 



APPENDIX A: DIFFERENTIAL FORMS: BRIEF 
REVIEW OF DEFINITIONS AND MAIN 
RESULTS 

A differential form w in a connected subset A of M 2 is 
formally written as: 

cj = M.(xi, x 2 )dxi + J\f(xi, X2)dx2 (Al) 

where M.{x\,x 2 ),N{x\,x 2 ) are two functions on A and 
(x\,X2) are the coordinates in IR 2 . (35j 
Given a curve tp : [sq, Si] — ► A, 

*Ks) = (M«),Ms)) (A2) 
the integral of lu along the curve -0, is defined as: 



u= / [M(mWi(t)+M(mW2(t)}dt (A3) 



where ip'^ 2 are the derivatives of ipi s2 . 

A differential form lu is said to be exact, if there exist 
a function G : A — > K such that: 



LU = dG 



(A4) 



that is 
dG 



dG 



(xi,X 2 ) = M(X1,X 2 ), g^-(xi,X 2 ) = M{X\,X2) 

(A5) 
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G is called a primitive for the differential form. 



Similarly we also get: 



Let H(A,B) be the set of all curves connecting the 
point A = (ai,a 2 ) to the point B = (61,62) m A. A 
differential form is exact if and only if for any couple of 
points A and B in A and curves ip and <p> in S(A, B), it 
is 



(A6) 



The integral of an exact differential form along any 
curve ip connecting A to B does not depend on the curve 
ip, but only on the ending points, and is given by: 



uj= dG = G(B) - G(A). (A7) 



The following statement also holds: A differential form 
is exact if and only if its integral along any closed curve 
is null. 

If the functions M. and M are of class C 1 (i.e. they are 
differentiable), then a necessary condition for the form w 
to be exact is that: 



dM dN 



(A8) 



3x2 dxi 

In this case the differential form lu is said to be closed. 

APPENDIX B: PROOF OF THEOREM 1 

It is known [l5[ that in one-dimensional systems con- 
fined in a [/-shaped potential, the period r of the orbit 
is equal to the derivative of the area $ in phase space 
enclosed by the orbit with respect to energy. 



<9$ 

tm 



(Bl) 



One simple way to prove this relation is by expressing 
the area as <&{E,V) — J dpdx9(E - H{x,p;V)), where 
9(x) is Heaviside step function [0(x) = 1 if x > 0, 9(x) = 
if x < 0]. Taking the derivative with respect to E, 
and using the relation S(x) = d6(x)/dx, gives r (see Eq. 
([27)). Using Eq. JbTJ) and Eq. (19]) we obtain: 



dS , t 
dE =kB *- 



From Eq. (f23|) , we obtain the relation: 



2*0* = - 

T 



from which, using Eq. (|B2|) . we get: 

dS_ k B 

OE ~ 2(K) t 



(B2) 



(B3) 



(B4) 



55 _ k B I dip \ 
Using Eqs. (|13I14|) we obtain: 



dS_ 
~dE 

&s_ 

dV 



I 

T 
P 
T 



(B5) 

(B6) 
(B7) 



APPENDIX C: PROOF OF THEOREM 2 

The proof of Theorem 2 makes use of the multidimen- 
sional version of Eq. (|BI[) . that is: 



O 



9$ 



N 



N 



dE 



(CI) 



This can be proved, in a similar way, by expressing <&at as 
J dqdp6(E — H(q,p;V)) and using the relation 5(x) = 
d0{x) / dx. The equipartition theorem [l7j | 



2(A') A 
3A 



(C2) 



is the generalization of Eq. (|B3[) to many dimensions. 
Using (J22J) and (JUT) with Eq. J57]) we get: 



95 



iV 



3Afc B 

In a similar way we also get: 

ds N 



dV 



3Nk B /d<p N ' 
2(A W ) AI \ dV j \ 



(C3) 



(C4) 



Using Eqs. (|34l35p with the ergodic hypothesis, we 
finally arrive at: 



(C5) 
(C6) 



dS N 


1 


dE 


Z Tn 


dS N 


Pn 


dV 





APPENDIX D: PROOF OF ADIABATIC 
INVARIANCE OF $ N 

We consider the time-dependent Hamiltonian 

H N (q, p; V(t)) = K(p) + ^(q; V(t)). (Dl) 



To prove that in Eq. ([38]) is an adiabatic invariant 
we first take the total time derivative of the Hamiltonian 
H N (q,p;V(t)) inEq. 02), 



dH N {q,p;V) dH N {q,p;V) dV 



dt 



dV 



dt 



(D2) 
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where the terms involving q and p cancel by Hamilton's 
equations. The derivative dV/dt changes slowly in 
time, but dHx/dt and dH^/dV can change rapidly be- 
cause of their dependence on q(t) and p(t). To eliminate 
the fast variables q, p we take the average of Eq. (|D2|) 
with respect to the microcanonical ensemble, which gives 



dH 



N 



dt 



dH 



N 



dV 



dV 
~dt' 



(D3) 



APPENDIX E: PROBLEMS 



Consider the following Hamiltonian of a one dimen- 
sional harmonic oscillator with angular frequency V (see 
Fig. [J) 



By Liouville's theorem [36| the average on the left-hand 
side of Eq. (|D3|) is 



dH 



N 



dt 



dE 
~dt' 



(D4) 



The microcanonical average in Eq. (|33p on the right- 
hand side of Eq. l[D3]l is 



dH 



N 



dV 



J ■ ■ J cMp p N (q, p, E, V) 



dH N (q,p,V) 



dV 



i a$ 



N 



n N dv 



(D5) 



where $at and SIn are given in Eqs. (|32|) and |38|) . re- 
spectively. Substituting Eqs. (jD4| and (|D5[) into Eq. 
(ID3|) and using £1^ = <9$ n / dE, we obtain 



d<P N dE d$ N dV 



dt 



dE dt dV dt 



= 0, 



(D6) 



which shows that <&jv is constant and therefore an adia- 
batic invariant. 



H(x,p; V) = 



2m 



(El) 



(a) Calculate the area ${E, V) enclosed by the trajectory 
of energy E and angular frequency V . Using Eq. (|B1|) 
check that the period of the orbit is, as expected, given 
by t(E,V) = 2tt/V. 

(b) Using Eqs. QH HU) show that k B T(E, V) = E, 
P(E,V) = -E/V 

(c) Show that the differential form dE + PdV, with 
P(E, V) as in (b) is not exact. (Hint: show that 
Eq. (|A8[) i s n °t satisfied.) Show that the integral 
of dE + PdV over the rectangular path with corners 
(E Q , V ), (E , Vi), (Ei,Vi), (E^Vo), and E ^ E x V + 
V±, is not null. 

(d) Consider the differential form ui = (l/T)dE + 
(P/T)dV, with P(E, V), T(E, V) as in (b). Find a prim- 
itive function S(E, V) for ui. Show that, apart from an 
additive constant, it is S(E, V) — log$(£', V), as dic- 
tated by Theorem 1. Check that Eq. (|A8[) is satisfied. 
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